Electron relaxation in a double quantum dot through two-phonon processes 
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We theoretically study the relaxation of electron orbital states of a double quantum dot system 
due to two-phonon processes. In particular, we calculate how the relaxation rates depend on the 
separation distance between the quantum dots, the strength of quantum dot confinement, and the 
lattice temperature. Enhancement of the rates by specific inter-dot distances and lattice tempera- 
tures, and the relative strength of different scattering channels are discussed. Our results show that 
although at low temperatures (T ~ 1 K) two-phonon processes are almost four orders of magnitude 
weaker compared to one-phonon processes in relaxing electron orbital states, at room temperature 
they are as important as one-phonon processes. 

PACS numbers: 73.21. La, 71.38.-k, 03.67.Lx, 85.35.Be 
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I. INTRODUCTION 

Semiconductor and superconductor nanostructure 
based quantum computing has attracted wide spread 
attention in recent years. ^ Specifically, single electrons 
trapped in a semiconductor double quantum dot (QD) 
have been suggested as candidates of quantum bits 
(qubits), where the location of the electron in the two 
dots represents the two states of the qubit iSiSiiiSiLi The 
proposal of these so-called charge qubits has prompted 
careful experimental investigations^iS during the past two 
years, which in turn motivated several detailed theoreti- 
cal studies of the decoherence properties of these charge 
qubit SiifliiLiSii? So far all the calculations on electron re- 
laxation due to electron-phonon interaction have focused 
on one-phonon processes, specifically single-phonon emis- 
sion process. Here we calculate the relaxation rates due 
to two-phonon processes and investigate whether there 
exist any regime where two-phonon processes might be 
as important as one-phonon processes to charge qubit 
decoherence. 

Two-phonon processes in a single three dimensional 
isotropic GaAs QD were studied in the context of electron 
relaxation in photoluminescence experiments more than 
a decade ago,-'^'* where processes of the most interest in- 
volve a longitudinal acoustic (LA) and a longitudinal op- 
tical (LO) phonon in the bulk phonon approximation.^^ 
Here we study the two-phonon processes in the context of 
relaxation of double dot charge qubits, where the qubit 
energy splitting (< 1 meV) is much smaller than that 
studied before (> 36 meV, when LO phonons are gener- 
ally involved) so that our focus will be on two-acoustic- 
phonon processes. In these processes phonon density of 
state, qubit energy splitting, and lattice temperature to- 
gether lead to interesting behaviors in the overall relax- 
ation rates and relative strength of the different scatter- 
ing channels. 

The paper is organized as follows. In Section II we 
first briefly describe the electronic states we study and 
the electron-phonon coupling in our system. We then 
give the expressions for the two-phonon relaxation rates. 
In Section HI we present our calculation results, showing 



how the relaxation rates depend on double dot param- 
eters such as inter-dot distance and single dot confine- 
ment, and how these rates vary with temperature. We 
also discuss the physical mechanisms behind these ob- 
tained behaviors. In Section IV we draw some conclu- 



II. THEORETICAL FORMALISM 

Our model system consists of two coupled QDs sep- 
arated by a distance of 2a, each of which is described 
by a two-dimensional (2D) harmonic wellii^ The single- 
dot one-electron wavefunctions are 2D harmonic oscilla- 
tor functionsi^ii^ and are described in terms of the prin- 
cipal quantum number n = 0,1,2,... and the angular 
momentum quantum number m — 0,±1,±2, ... as 
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where (p^) are the Laguerre polynomials, and p = 

|r|||// is a scaled radius, with I = y^fi/m*u!o. The corre- 
sponding eigenvalues are 
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Along the growth direction we assume an infinite quan- 
tum well (QW) confinement, so that the electron wave- 
function has the form 
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For two QDs that are horizontally coupled, we use a sim- 
ple in-plane confinement of two parabolic wells separated 
by an inter-dot distance 2a: 

Vc ='^m*u}l imn{{x - af + , {x + af + y'^} . (4) 

The single electron wavefunction for the lateral direc- 
tion can in general be expressed as a superposition of the 
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single-dot wavefunctions: 



(5) 



and the total wavefunction of the system of the coupled 
QDs is 



*(r) = (r||)Vz(z) 
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In the present study, the in-plane wavefunctions for the 
coupled-QD are calculated numerically by direct diago- 
nalization, using reasonable parameters of a GaAs QW. 

Electrons in GaAs interact with both acoustic and op- 
tical phonons. However, here we do not consider contri- 
butions from the optical phonons to the electron relax- 
ation due to the small energy splitting in the double dot 
system we study. In this work, we calculate two-phonon 
relaxation rates caused by both deformation potential 
and piezoelectric interactions. The Hamiltonian that de- 
scribes these interactions is given by: 
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Here emission of two phonons (-I-LA-I-LA) , absorption of 
two phonons (-LA-LA), and emission and absorption of 
one phonon each (-f LA-LA, including both overall emis- 
sion and absorption) are indicated by subscripts , 

and -I — . The later two processes vanish at T = K. 
However, since experimental temperature is often in the 
same order as the qubit energy splitting, we include all 
two-phonon processes in our calculation. Indices i and 
/ represent the initial and final electronic states, which 
are the ground and first excited double dot states respec- 
tively. Index s refers to intermediate electronic states 
which are among the 13 higher energy states we include 
in the present calculation (the summation over s excludes 
the initial and final states, and the sum has converged). 
The matrix elements are calculated as in Ref. Notice 
here that both deformation potential and piezoelectric in- 
teractions are included in the matrix element — they can- 
not be treated separately as in one-phonon processes. iVk 
is the Bose distribution function for the k phonon mode 
with energy _Ek — ?iwk- The integrals over k and q in 
Eq. lfTn|) to (fT^ are calculated by Monte Carlo code. 



III. RESULTS AND DISCUSSIONS 



where /9„ is the mass density of the host material, ujq is 
the frequency of the phonon mode with wave vector q, 
V is the volume of the sample, Oq and aLq are phonon 
annihilation and creation operators, and p(q) is the elec- 
tron density operator. The interaction strength A^(q) is 
defined by 



Miq) = D\q\ + iMx{q), 



(8) 



where the first term represents the deformation poten- 
tial interaction with deformation constant D, and the 
second term, which is imaginary, describes the piezoelec- 
tric interaction. For zincblende crystals (e.g. GaAs), the 
piezoelectric term Ai\{q) takes the form 

-^r(q) = 2e ei4 {qxqy£,z + qyqzCx + q-xqz^y) (9) 

Using second-order perturbation theory, scattering 
rates due to the emission and/or absorption of two LA 
phonons can be obtained 
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At low temperatures, which are the operating temper- 
ature of the current generation charge-qubit-based ar- 
chitectures, two-phonon absorption processes should be 
much weaker than emission processes. Therefore we first 
evaluate the rates for an electron in the excited qubit 
state (the first excited state of the double dot as the ini- 
tial state) to relax to the ground qubit state (final state) 
via two-phonon processes. Two channels are included 
here: (a) emission of two LA phonons (+LA+LA) and 
(b) emission of a high energy LA phonon and absorption 
of a lower energy LA phonon (-I-LA-LA). Throughout the 
paper the QW width is fixed at 2Lz = 10 nm and the 
material parameters are taken from Bruus et al.^^ 

In Figs, n and 121 we present the calculated relaxation 
rates due to the two overall emitting processes as func- 
tions of half inter-dot distance a and single dot confine- 
ment energy. For both figures the lattice temperature 
is fixed at T = 1 K, a usual operating temperature for 
semiconductor charge qubits. The results presented in 
Fig. n] have some interesting features. The most promi- 
nent is that the order of magnitude for both processes is 
much smaller compared to the one-phonon processes^'^. 
Therefore at low temperatures it should be sufficient to 
consider only the single-phonon processes when studying 
phonon-induced charge qubit decoherence. 

Both curves in Fig. ^ display a clear peak at an in- 
termediate inter-dot distance. This feature is an in- 
terplay between phonon density of state and the mag- 
nitude of the electron-phonon matrix element. When 
inter-dot distance increases, the qubit energy splitting 
decreases monotonically. This leads to a decreasing 
phonon density of state as it varies as as a function of 
phonon frequency. On the other hand, the magnitude of 
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FIG. 1: Relaxation rates of an electron in a double dot due to 
two different two-phonon processes (+LA+LA and +LA-LA) 
as functions of the half interdot distance a. The confinement 
strength is hui — 3 meV, lattice temperature T = 1 K, and 
QW width 2Lz = 10 nm. The oscillatory behavior of the 
+LA+LA process at larger a is an artifact of the integration 
process. 
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FIG. 2: Relaxation rates of an electron due to two different 
two-phonon (overall emission) processes (-fLA+LA and -l-LA- 
LA) as a function of the confinement strength. The rates 
versus the energy splitting Ae are shown in the inset. The half 
interdot distance is a = 25 nm, and the lattice temperature 
is T = 1 K. 



electron-phonon matrix element increases with increasing 
ai^ These two opposite trends dictate that a maximum 
should develop at an intermediate inter-dot distance for 
the electron relaxation rate. 

One additional feature of Fig. ^ is the crossover be- 
tween the two curves. The physics behind this crossover 
is even more complicated than the maxima of the two 
curves at an intermediate value of a. Qualitatively, 
for +LA+LA processes, both phonons involved have to 
have energies smaller than the qubit energy spUtting, 
while for -f LA-LA processes the emitted phonon has to 
have energy larger than the qubit energy spUtting. As 
bulk acoustic phonons have a density of state propor- 
tional to Lj^, the +LA-LA processes should be favored 
over +LA+LA processes. On the other hand, -fLA- 
LA processes involve absorption of one phonon, which 
is thermally constrained so that its energy has to be 
relatively small (thus with a relatively small phonon 
density of state) compared to the thermal energy (in 
the current case ^ 0.1 meV). This argument favors the 
-l-LA+LA processes. Furthermore, we mentioned above 
that the electron-phonon interaction matrix element in- 
creases with increasing a. Therefore, which process 
yields a higher relaxation rate is determined by consid- 
ering all three factors of lattice temperature, electron- 
phonon matrix element, and phonon density of state. 

The results presented in Fig. [21 are obtained at a fixed 
interdot distance of 2a — 50 nm. Notice that here the 
maxima of the two relaxation rate curves are mostly 
determined by the qubit energy splitting (as shown in 
the inset) and the phonon density of state consideration. 
The matrix element plays a less important role because 
it is not as sensitive to the single-dot confinement en- 



ergy as it is to inter-dot distance. Again the two rates 
have a crossover, with -t-LA-l-LA processes faster at larger 
qubit energy splittings while -f-LA-LA processes faster for 
smaller qubit energy splittings. The physical mechanism 
is similar to that in Fig. ^ 

In Fig. O we present the electron relaxation rates as 
functions of the lattice temperature for all the one- 
phonon and two-phonon processes. At T = 1 K, the 
strongest two two-phonon processes are the -|-LA-|-LA 
and -l-LA-LA overall emission process, but the corre- 
sponding relaxation rates are more than three orders of 
magnitude smaller than the one-phonon emission pro- 
cess. The other two two-phonon processes are thermally 
suppressed because of the involvement of absorption of 
high energy phonon(s). On the other hand, at room tem- 
perature (corresponding to a thermal energy in the order 
of 30 nieV, much larger than the qubit energy splitting of 
< 1 meV), the strongest two-phonon process, the -l-LA- 
LA (overall emission) process is only one order of magni- 
tude smaller than the one-phonon processes (both emis- 
sion and absorption). In other words, if the charge qubit 
states of a double quantum dot is to be used for some 
form of quantum information processing at high temper- 
atures such as room temperature, two-phonon relaxation 
processes should also be taken into consideration. 

At intermediate temperatures the relaxation rates 
have crossovers between the one-phonon-emission-one- 
phonon- absorption processes and the two-phonon ab- 
sorption or emission processes. The physical mecha- 
nism is similar to what we have discussed for Figs. Q] 
and 12 in that the thermal suppression of the higher en- 
ergy phonons (which are only involved in +LA-LA and 
-LA-I-LA processes and have higher density of state) is 
gradually removed as temperature increases. At room 
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FIG. 3: Relaxation rates of an electron due to all the two- 
phonon processes [-I-LA-I-LA, +LA-LA (overall emission), - 
LA-LA (two-phonon absorption), and -LA+LA (overall ab- 
sorption)] and one-phonon processes [-I-LA (one-phonon emis- 
sion) and -LA (one-phonon absorption)] as functions of the 
lattice temperature. Here the confinement strength is ?itj = 3 
meV, and the half interdot distance is a = 25 meV. 



temperature the -|-LA-f LA and -LA-LA processes have 
almost the same relaxation rates because the thermal fac- 
tors in Eqs. (|10|l and are approximately the same for 
phonons with energy smaller than 1 meV (recall that for 
two-phonon emission or absorption processes each of the 
phonons involved should have energy smaller than qubit 



energy splitting). Similarly, -t-LA-LA (overall emission) 
and -LA-fLA (overall absorption) processes also have 
similar rates at high temperatures. The discrepancy is 
mostly due to the involvement of higher energy phonons 
in these processes, for which the factors Nk and {Nk + 1) 
could be sufficiently different. 



IV. CONCLUSIONS 

In this study we have investigated two-acoustic-phonon 
processes (including both phonon emissions and absorp- 
tions) induced charge qubit relaxation in a semiconductor 
double quantum dot. Our results show interesting depen- 
dence of the relaxation rates on the lattice temperature 
and system configuration parameters such as inter-dot 
distance and single dot confinement energy. We have 
found that acoustic phonon density of state, electron- 
phonon coupling matrix element, and phonon thermal 
distributions together lead to crossovers between differ- 
ent two-phonon scattering channels. In the context of 
charge qubit based quantum information processing, the 
two-phonon processes are much weaker than the one- 
phonon processes at the relevant low temperatures, al- 
though their magnitudes are similar at the room temper- 
ature. 
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